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Using the WMAP cosmic microwave background data it is shown that collisions between Alfven 
wave packets and the cascades generated by these collisions (the Iroshnikov model) can determine 
the photon temperature fluctuations for arcminute scales on the last scattering surface. 
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Cosmic electromagnetic fields are expected to be gen- 
erated by the cosmological phase transitions (electroweak 
and QCD) in a wide interval of scales before 
the recombination time. At earliest times, the magnetic 
fields are generated by particle physics processes, with 
length scales typical of particle physics. It is shown in Q 
that turbulence with its cascade processes is operative, 
and hence the scale of magnetic fields is considerably 
larger than would be the case if turbulence were ignored. 
In particular, it is shown in |3j,[4|,|g] that rotational 
velocity perturbations, induced by a tangled magnetic 
field can produce significant angular scale anisotropies in 
cosmic microwave background (CMB) radiation through 
the Doppler effect. The conclusions are relevant to 
arcminute scales 0. For very large cosmological scales 
(larger than 10°), on the other hand, the Anderson lo- 
calization can make it impossible for the electromagnetic 
fields to propagate [l^. Magnetohydrodynamic (MHD) 
turbulence is characterized by a competition of two 
processes, Alfven wave packets collisions and swirling 
motions. It is known |4| that on the last scattering 
surface the nonlinear Alfven wave mode survives photon 
(Silk) damping on the arcminute scales, while the 
general swirling motions (as well as the compressional 
modes) are effectively dissipated on these scales (it is 
significant for further consideration that the Alfven wave 
modes induce specific rotational velocity perturbations 
on the last scattering surface 0,0], 0,133). Therefore, 
the cosmic baryon-photon fluid becomes dominated by 
Alfven waves on the arcminute scales just before the 
recombination time. In the present paper we will show 
that specific statistical properties of such Alfven wave 
dominated fluctuations are consistent with the new 
(WMAP) arcminute CMB data. 

The incompressible magnetohydrodynamic equations 
can be written in terms of the Elsasser variables 

z ± =v±B (1) 

as 

a t z ± +z =F -Vz ± = -VP + ^+Az ± +^_Az T , (2) 
where v+ = \(y + n), V- = -^{y — rj), P is the total pres- 



sure, v and r\ are coefficients of hydro and magnetic diffu- 
sion respectively. The equations are given in convenient 
nondimcnsional form using units B / Bq — ► B, u/uq —> u, 
Uq = BqK^ttp) 1 / 2 and B is a typical magnetic field in- 
tensity. 

Scaling of the structure functions of the Elsesser vari- 
ables 
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is used as an effective tool to study their dynamics. 

A first attempt to describe magnetic turbulence dom- 
inated by Alfven waves was made in [l9| . In the in- 
compressible fluid, any magnetic perturbation propa- 
gates along the magnetic field line. Since wave pack- 
ets are moving along the magnetic field line, there are 
two possible directions for propagation. If all the wave 
packets are moving in one direction, then they are sta- 
ble. Therefore, the energy cascade occurs only when the 
opposite-travelling wave packets collide, and only colli- 
sions between similar size packets are taken into account 
in the Iroshnikov model. The following amount of en- 
ergy: AE ~ (vf /1)(1/Va) is released at collision of the 
two wave packets of the same size I. The energy change 
per collision and the duration of the collision are respec- 
tively vf (vi/Va) and At ~ I/Va, where Va is the Alfven 
speed. Total number of collisions for the cascade can be 
estimated as v 2 j AE. Hence, the energy cascade time t; 
is 



n ~ {v 2 /AE) 2 At 
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That is the cascade time is (Va/vi) times longer than the 
eddy turnover time (l/v{). The constancy of energy cas- 
cade:, (vf)/(lVA) = constant, is assumed in the model, 
which results in 

vf ~ I. (5) 
It then follows from (4) and (5) 

n ~ I 1 ' 2 . (6) 

There exists a general representation for Q p (3) |2(| 
C P = p/g(l - x) + C [l - (1 - x/Co) p/s ], 
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where g is related to the basic scaling Svi ~ I 1 ' ' 9 (in the 
Iroshnikov model (5) gives g = 4), x is the scaling of 
dynamic time scale of the most intermittent structures, 
T\ ~ P (in the Iroshnikov model (6) gives x = 1/2), Co 
is the co-dimension of these structures with the spatial 
dimension d, Cq = 3 — d (in the Iroshnikov model these 
structures are known to be micro-sheets, i.e. d — 2 and, 
consequently, Cp = 1). Then, for the Alfven wave domi- 
nated model |21j 



Cp = | + l-(l/2) p/4 . 



(7) 



It is very significant that for this model (4 = 1 (cf. basic 
scaling (5)). Using this fact, the so-called extended self- 
similarity (ESS) can be introduced in the form J2jj 

flz^x + r) - z±(x)n ~ <|z ± (x + r) - z ± (x)| 4 )^. (8) 

The remarkable property of the ESS is that the multi- 
scaling (8) can survive even if the original multiscaling 
(3) does not exist. Due to the local (small-scale) isotropy 
of the fluctuations in the model we suppose that the ESS 
has the same £ p for any vector component of the space 
difference: z (x + r) — z ± (x). 

We use the ESS and the model (7) to check whether the 
cosmic microwave background (CMB) data, obtained re- 
cently by the WMAP space mission, support the Alfven 
wave domination on the arcminute scales (cf Introduc- 
tion). For this purpose we calculated moments (AT,?), 
for the space differences of the CMB temperature fluctu- 
ations AT r = |T(x + r) — T(x)|, and then checked the 
multiscaling 



(AT r p ) ~ (AT r 4 ) c " 



(9) 



(cf. (8)). 



The motion of the scatterers imprints a temperature 
fluctuation, 5T, on the CMB through the Doppler effect 
@,H (see also and 0) 



T 



J g(L)n ■ v b (x) dL 



where n is the direction (the unit vector) on the sky, v b 
is the velocity field of the baryons evaluated along the 
line of sight, x = Ln, and g is the so-called visibility. 
It should be noted that, in a potential flow, waves per- 
pendicular to the line of sight lack a velocity component 
parallel to the line of sight; consequently, generally there 
is no Doppler effect for the potential flows. The same is 
not true for vortical flows, since the waves that run per- 
pendicular to the line of sight have velocities parallel to 
the line of site [23 ( see 0j0 about the rotational veloc- 
ity perturbations induced by the Alfven-wave modes on 
the last scattering surface). One can look at the modula- 
tion from another (though similar) point of view [6|. In 
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FIG. 1: a: Logarithm of moments of different orders (|AT,-| P ) 
against logarithm of (|AT r | 4 ) for the cleaned and Wiener- 
filtered WMAP1 data set. The straight lines (the best fit) 
are drawn to indicate the scaling (9); b: the same as in figure 
la but for WMAP2 data set. 



general, vector perturbations of the metric have the form 







Bi 



Bj Hi j -)- H 
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where B and H are divergence-free, 3D vector fields sup- 
posed to vanish at infinity. The authors of 6] introduced 
two gauge invariant quantities 



a = H B and 



v-B 



which represent the vector contribution to the perturba- 
tion of the extrinsic curvature and the vorticity (cf. (6) 
and (13)). The general form of the CMB temperature 
fluctuations produced by vector perturbations is 



8T 



(vec) 
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where Z = z — a is a gauge-invariant generalization 
of the velocity field, and the subscripts dec and 
denote the decoupling epoch (zdec 3> 1100) and today 
respectively. We see from this equation that, besides 
the Doppler effect, Alfven waves gives rise to an inte- 
grated Sachs- Wolfe term. However, since the geometric 
perturbation a decays with time, the integrated term is 
dominated by its lower boundary and just cancels a in 

We used the WMAP data cleaned from foreground 
contamination and Wiener filtered |23j (the original tem- 
perature is given in /iK units). Wiener filtering sup- 
presses the noisiest modes in a map and shows the signal 
that is statistically significant. However, even after the 
cleaning and filtering some data points with the largest 
magnitudes seem to be suspicious. Therefore, we also 
made several magnitude-threshold cutoffs to check the 
stability of our calculations to these cutoffs. 
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FIG. 2: The scaling exponents £ p , corresponding to the 
WMAP data (circles - WMAP1, triangles - WMAP2). The 
dotted straight line corresponds to the Gaussian distributions 
and the dashed curve corresponds to the Alfven-wave domi- 
nated model calculations (7). 



FIG. 3: a: Logarithm of the temperature gradient moments 
(X?) calculated for the WMAP1 data set against Ins. The 
straight lines (the best fit) are drawn to indicate the scaling 
in the log-log scales; b: the same as in figure 3a but for 
WMAP2 data set. 



Figure la shows the results of the calculations in the 
form suitable for the multiscaling of the type (9). The 
straight lines in this figure (the best fits) correspond to 
the multiscaling (9). The cutoff in this case excluded 
about 1% of the data points (we will call this data set 
WMAP1). Figure lb shows analogous data with a cutoff 
which excluded about 10% of the data points (we will 
call this data set WMAP2). 

Figure 2 shows the exponents ( p (open circles corre- 
spond to WMAP1 data set and triangles to WMAP2) 
extracted from figures la,b as slopes of the straight lines. 
We also show in figure 2 as dotted line the dependence 
of f p on p for the Gaussian distributions (( p — p/A for 
any Gaussian distribution), and we show the model de- 
pendence (7) as dashed line. 

One can see that starting from p ~ 6 the data depart 
systematically from the Gaussian straight line and 
follow quite closely to the model curve (indicated by the 
dashed line) predicted for the Alfven-wave dominated 
turbulence (7). 

To give an additional support to the correspondence 
just noted between the WMAP data and the Alfven wave 
dominated model, let us introduce a gradient measure for 
the cosmic microwave radiation temperature T fluctua- 
tions 



Xr 



L„(vT) 2 dv 



(10), 



where v r is a subvolumc with space-scale r. Scaling laws 
of this measure, such as 



(X P r) 
(Xr)? 



(11) 



are an important characteristic of the temperature dissi- 
pation rate [24[ ■ The exponents (j, p can be related [24| to 
the exponents ( p by the equation 



A'p — 1 dp 



(12) 



that allows us to check the model equation (7) also 
through the scaling (11) of the moments of the gradient 
measure. 

Technically, using the cosmic microwave pixel data 
map, we will calculate the cosmic microwave radiation 
temperature gradient measure using summation over 
pixel sets instead of integration over subvolumes v r . The 
multiscaling of type (11) (if exists) will be then written 
as 



(X P s) 
(Xs) P 



(13) 



where the metric scale r is replaced by number of the 
pixels, s, characterizing the size of the summation set. 
The Xs is a surrogate of the real 3D dissipation rate Xr- 
It is believed that the surrogates can reproduce quanti- 
tative multiscaling properties of the dissipation rate [24| . 
Since in our calculations (x s ) is independent of s, we will 
calculate the exponents fi p directly from the scaling of 

Figure 3a shows scaling of the CMB temperature gradi- 
ent moments (xf) calculated for the WMAP1 map. The 
straight lines (the best fit) are drawn to indicate the scal- 
ing in log-log scales. Figure 3b shows the results of anal- 
ogous calculations produced for WMAP2 data set. 

Figure 4 shows the exponents \i p extracted from figure 
3a (circles) and from figure 3b (triangles). Dashed line 
in figure 4 corresponds to the model calculations (eqs. 
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FIG. 4: The scaling exponents u p corresponding to the 
WMAP1 data set (circles) and for the WMAP2 data set (tri- 
angles). The dashed curve corresponds to the Alfven -wave 
dominated model calculations (7), (12). 



the non-Gaussianity: using the fluctuations of the tem- 
perature (not their space increments, (9)) and the fact 
that non-Gaussianity does not make its appearance in 
the ESS up to the 6 moment even if space increments are 
considered. 

Finally, it should be noted that we cannot exclude 
other contributions to high-order structure functions 
which may be produced by non-linear effects due to 
gravity (e.g. the Rees-Sciama effect), but the good 
quantitative correspondence to turbulence (both for 
the structure functions and for the gradient measure), 
observed in the figures 2 and 4, indicates that the 
turbulence may be a dominating factor here. 

The authors are grateful to C.H. Gibson for discus- 
sions, and to Tegmark's group and to the NASA Goddard 
Space Flight Center for providing the data. 



(7),(12)). 

The results presented in the paper can be considered 
as a tentative indication of the existence of the 

considerable magnetic fields at the recombination time. 
The extended self-similarity (figure 2) shows clear non- 
Gaussian character of the high moments which, together 
with the closeness of the low moments to the Gaussian 
behavior, can shad light on the longstanding discussion 
about Gaussianity of the small-scale cosmological fluctu- 
ations. The theoretically predicted survival of the Alfven 
waves on the last scattering surface at the arcminute 
scales and, moreover, domination of their interactions 
(the cascade collisions of the Alfven wave packets) on 
these scales can be also considered now rather plausi- 
ble. WMAP sky maps produced by the WMAP team 
were analyzed by the team itself. They found the CMB 
temperature fluctuations (i.e. ST — T — (T)) to obey 
Gaussian statistics :25|. Analogous situation occurs in 
fluid turbulence (near-) Gaussianity of velocity fluctua- 
tions and non-Gaussian multiscaling of the correspond- 
ing space increments |2^,|2^. The (near-) Gaussianity of 
a random field fluctuations themselves does not contra- 
dict to a pronounced non-Gaussian multiscaling of cor- 
responding increments of the field, which can be seen 
for sufficiently high order moments. This takes place 
in classic fluid turbulence for velocity field and, appar- 
ently, for the CMB temperature fluctuations (probably 
just due to turbulence modulation). In particular, it was 
shown for fluid turbulence (2(| , [27| that knowledge about 
the spectrum for the Gaussian processes is certainly in- 
sufficient to reproduce the observed multiscaling of the 
increments. We suppose that there are two reasons why 
the previous studies of the CMB maps failed to detect 
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